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A theory of the local magnetic response of a nonmagnetic impurity in a doped antiferromagnet, as relevant to 
the normal state in cuprates, is presented. It is based on the assumption of the overdamped collective mode in the 
bulk system and on the evidence, that equal-time spin correlations are only weakly renormalized in the vicinity of 
the impurity. The theory relates the Kondo-like behavior of the local susceptibility to the anomalous temperature 
dependence of the bulk magnetic susceptibility, where the observed increase of the Kondo temperature with 
doping reflects the crossover to the Fermi liquid regime and the spatial distribution of the magnetization is given 
by bulk antiferromagnetic correlations. 
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One of the open theoretical questions in cuprates is the un- 
derstanding of a well established experimental fact that non- 
magnetic impurities have a strong effect on superconduct- 
ing as well as on the normal-state properties of cuprates ll|]. 
Prominent example is the substitution of Cu in CUO2 planes 
by Zn-^+ which acts as a spinless impurity. In contrast to the 
naive picture of a weak scatterer, Zn-^+ impurity depresses 
superconductivity already at small concentration and repre- 
sents a strong scatterer for transport properties. In this con- 
tribution we mainly address the magnetic effects of an im- 
purity in the non-superconducting phase. NMR experiments 
show that Zn induces local magnetic moments on the near- 
est neighbor (n.n.) copper sites |2], as well as on more dis- 
tant copper neighbors |3l |4|]. The observed local suscepti- 
bility Isl la] is well accounted for by the Curie-Weiss form 
Xioc oc l/(r + 8), whereby 8 is independent of the impu- 
rity concentration but reveals a clear dependence on doping. 
Namely, in the underdoped YBaCuaOe+a; (YBCO) the behav- 
ior is nearly Curie-like, with 8 '^ 0, whereas in the overdoped 
YBCO 8 shows a strong increase with doping |7]. In anal- 
ogy with the Kondo effect in metals 8 has been interpreted 
as a relevant Kondo temperature. Analogous effects in doped 
cuprates have been established for a Li+^ impurity which is 
also nonmagnetic |7]. Since it has a different valence the sim- 
ilarity indicates that the additional hole is not trapped near 
the impurity. By a multi-nuclei NMR imaging it has become 
also increasingly clear that the spatial distribution of the mag- 
netization around the impurity reflects the antiferromagnetic 
(AFM) correlations of the bulk system J3|]. 

From the point of theoretical description, it seems well es- 
tablished that a nonmagnetic impurity as, e.g., Zn^+ can be 
incorporated into a microscopic electronic model of CUO2 
planes in cuprates by introducing, e.g., an inert empty site 
into the t-J model relevant to cuprates fST, or a local site 
with a very different local energy within the Hubbard model 
|9, 10]. Impurity-induced moments and a Curie-like suscep- 
tibility have been established in disordered spin systems, e.g., 
spin ladders ] 11], in the presence of the gap in the spin ex- 
citation spectrum. An analogous treatment, assuming a spin 
gap in underdoped cuprates, also leads to the appearance of 



a Curie-type susceptibility of an unpaired spin 11211 . Kondo 
effect based on the t-J model and on the spin-charge sepa- 
rated ground state has been used also to explain strong trans- 
port scattering on the impurity site ]13]. The 2D Hubbard 
model with impurities has been analysed using a renormal- 
ized random-phase-approximation (RPA) fy\ to describe the 
Knight-shift data in Zn and Li substituted YBCO. Numerical 
studies confirmed the existence of an induced Curie suscep- 
tibility at low doping ] 10], whereas variational Monte Carlo 
simulation established equal-time AFM correlations around 
the impurity 11411 . 

Nevertheless, the understanding of magnetic effects of such 
an impurity is far from satisfactory. It is clear that in the nor- 
mal state of cuprates there is at most a pseudogap in under- 
doped regime, whereas at higher doping spin excitations are 
gapless and the usual arguments for a Curie-type susceptibil- 
ity from unparied spins II2I 11311 become unfounded. There 
is also no explanation for the origin and the onset of finite 
and large Kondo temperature 8 > in the overdoped regime. 
Lacking is also the theoretical answer to the evident question 
whether the local magnetic response around the impurity is 
just the reflection of the bulk (as well anomalous) magnetic 
response, as evidenced by recent experiments |2|]. In any case, 
the T-dependence of magnetic response outside the spin gap 
regime has been so far treated only within the RPA ]9], which 
is generally not satisfactory enough to account for anomalous 
bulk properties. 

To address the above issues we present a generalization of 
the theory of spin response applied to explain anomalous bulk 
magnetic response ]15]. The novel input is the observation 
that equal-time spin correlations around the impurity are to a 
large extent unrenormalized. The consequence is that the local 
magnetic response is an image of the anomalous staggered 
susceptibility in a homogeneous system, which exhibits also 
the Curie-Weiss behavior Consequently, the onset of a finite 
Kondo scale 8 is related to the doping-driven crossover of the 
bulk spin dynamics from a non-Fermi-liquid spin dynamics to 
a Fermi-liquid one llfill . 

Let us start with an approach to the spin response of a ho- 
mogeneous doped AFM. The dynamical spin susceptibility 



Xq(w) = —{{S^; S^))ui can be expressed within the memory 
function formalism as 11511 



Xq(w) = 



-'7q 



UJ^ + UjAJqlLo) — 



(1) 



where7yq = -i([5iq,5^)]). 



'yq/XqandXq = Xq(tJ = 



0) is the static susceptibility, ujq represents the collective- 
mode frequency in the case of low damping 7q ^ M'' (wq) < 
Wq. However, in cuprates the collective mode is found to be 
overdamped throughout the normal phase, i.e., 7q > Wq. Fol- 
lowing the evidence from the analysis of microscopic models, 
such as the planar t-J model llsLllTll . we assume also con- 
stant rjq ^ rj and 7q ~ 7 in the region of interest near the 
AFM wavevector q ~ Q = (tt, tt). 77 is closely related to 
energy, so it is quite T-independent and only smoothly depen- 
dent on doping. In a doped AFM damping 7 emerges from the 
decay of a spin collective mode into electron-hole excitations 
and is also assumed to remain finite and large at low T — > 
in the normal state 1 17i]. Hence, the main T- variation enters 
Eq. O via Xq- 

The central idea of the theory of the anomalous scaling 
in doped AFM [15] is that a nontrivial dependence x<i{T) 
is driven by the fluctuation-dissipation relation for the equal- 
time correlations 
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There exists an extensive evidence from analytical 12011 and 
numerical work on the t-J model at finite hole concentration 
Ch> that Cq^Q -- C/(k2 + g^), q = q - Q, and the AFM 
correlation length ^ = 1/k saturates at low T. Similar conclu- 
sions emerge from an analysis of inelastic neutron scattering 
on YBCO lHH. The relation (|3 then leads to a strongly T 
dependent Xq^q (T) which is at heart of the non-Fermi liquid 
behavior observed in underdoped cuprates, in contrast to the 
usual Fermi liquid where Xq('^) is essentially T-independent. 
Let us first consider the behavior of Xq(r). Performing the 
high-T expansion in Eq. (|2jl we get 



1 r''^^ 

Cq-TXq+-y^ rf-^X»=TXq+^, 



(3) 



where we have taken into account the definition of 1]^ as the 
second frequency moment of the shape function x''{u!)/llj. 
The high-T expansion is thus consistent with the Curie- Weiss 
behavior xq = Cq/{T + 6q) where Oq = r7q/12Cq. 
Since in a do ped AFM, as represented, e.g., by the t-J model, 
Cq oc 1/ch 11611 we get quite small scale 0q ex Ch- On the 
other hand, the expansion is valid only for T > 0q. At T = 



we get XQiT = 0) = v/il^p) LiSl> where i 



76 



-2C and 



C = 7r7Cq/277. 

In Fig. 1 we present C'q/xq(T) as follows from the solu- 
tion of Eqs. ^M for fixed 77 = 0.6 i and 7 = 0.5 i 1 15] (note 
that t ^ 400 meV for cuprates) but varying k = 0.5 — 1.5 (in 
units of inverse lattice spacing). By way of Eq.(|3 k also de- 
termines k which enters the low-oj behavior of Xq{^)- For 



the range of k considered, k/k ~ 2 so that k roughly cor- 
responds to the experimental range of values in underdoped 
to overdoped YBCO ]18]. Note that the Curie-Weiss law 
is obeyed down to T ~ 6q. The deviation appears (to 
lower values for given parameters) on approaching T ^ 0. 
Qk ~ C'q/xq(T = 0) can be interpreted as the relevant 
Kondo temperature and its variation with k is presented in the 
inset of Fig. 1. On increasing k (note that k cx -^qT) we 
are facing quite an abrupt transition from a Curie behavior at 
T -^ to a Curie-Weiss variation with finite and rapidly in- 
creasing Qk > 1.16.1 . Such a behavior is also consistent with 
experimental results for l/xq(T) obtained from the analysis 
of NMR 1/T2g(T) relaxation in various cuprates lll6ll . 
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Figure 1: Cq/xq vs. T (both in units oft) for various k. The inset 
shows the variation of Kondo Qk with k. 

Let us turn to a doped AFM with an added nonmagnetic 
impurity, as relevant to Zn^+ replacing planar Cu in CUO2 
planes. To be more specific we have in mind the planar t- 
J model, where the impurity is represented as an empty site 
at the origin i = 0. In analogy with the homogeneous sys- 
tem, Eq.©, we consider the equal-time local correlations 
Cij — (SIS':) as an essential input. In general, C^ differ 
from correlations in an homogeneous system without any im- 
purity where C" ~ C°(Rj — Ri). (Furtheron we lable the 
homogeneous quantities by the superscript 0.) However, it 
appears characteristic for strongly correlated electrons in low- 
doped AFM that at least shorter-range Cy close to the im- 
purity deviate modestly from unperturbed C,° . This is partly 
plausible since, e.g., the n.n. correlations Cij are governed by 
the minimization of the exchange energy. On the other hand, 
insensitivity of Cij around the impurity can be interpreted as 
an effective spin-charge decoupling, where the empty site rep- 
resentsjust a free spinon not affecting nearby spin correlations 
112, 13]. 

In support of the above conjecture we perform an exact- 
diagonalization calculation of Cij within the t-J model with 
an empty site. As an example of a doped AFM we present 
in Fig. 2 results for a system A^ = 20 sites at T = with 



Nh = 2, 3 mobile holes and J/t = 0.3. We show some where we have for simplicity assumed the diagonal form of 
nonequivalent n.n. and next n.n. correlations Cij around the {Mij,r]ij) ~ X]a(*7-'^''M)(^a)*^a as well. In analogy with 
impurity. We notice that even on sites neighboring the im- the homogeneous system, the local magnetic response around 
purity at least shorter-range C^ are nearly the same as in the the impurity will be determined by the behavior around A ~ A 
bulk, or even enhanced as noted also by others L14I . for which Ca — max. In this region we assume constant 7a ~ 

7 and rjx ^ -q. Clearly, this is based on the assumption that 
damping is quite local and not affected significantly by the 
impurity. On the other hand, rjij can be explicitely calculated 
and is expressible in terms of equal-time correlations 11711 . 

In a system with a nonmagnetic impurity a spin polariza- 
tion around the impurity induced by a homogeneous magnetic 
field, gives rise to the susceptibility Xi on site i 
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Figure 2: Equal-time correlations Ctj around the impurity (denoted 
by ®) within the t-J model for J/t = 0.3 and dopings Ch ~ 2/20 
(upper) Ch = 3/20 (lower value) in units of 10^*. 

It is therefore plausible to treat a system with an impurity 
by assuming frozen correlations. I.e., we take that Cy = for 
either i = or j = 0, while Cij = C°, elsewhere. Note that 
at T = and in a homogeneous system on a bipartite lattice 
(with even N, Nh) the ground state is a singlet (S'tot — 0). In 
the case of a single impurity this leads to 
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-Ch, (4) 



since Y.J C% = (5f 5fot) = 0. Eq. gj reproduces for Ch < 1 
the proper moment of the impurity corresponding to S'tot — 
l/2L8,i2J. 

We further study the real-space local susceptibilities 
Xij(w) = —{{Sf; Sj))uj in analogy to the homogeneous case. 
The dynamical matrix xi^) can be generally expressed in 
terms of corresponding matrices in real space x, ry, M_ as 



X; 



"=E^*j =E^^(^5v)*«A, 



vx = 



E^ 



(8) 



In the case of a n.n. site Xi is directly related to the Knight 
shift on the ^Li impurity and of ^^Y near the Zn impurity |^. 
Another relevant quantity is the (average) uniform suscepti- 
bility X — X^A XaI^aP/-^ which yields the impurity-induced 
contribution Ax = X - x"=o- 

We first give an approximate solution to the impurity prob- 
lem via the perturbation calculation. Here, C_ = C^ + C' 
with the perturbatrive part Cqq = —1/4 and CL = C'q = 
— Cg.. The unperturbed eigenvectors are the homogeneous 
ones, Uq = exp{iqri)/\/N, and the lowest order calculation 
gives 



Av„ = 



1 

N 



E^ 






(9) 



where fi ~ ^1/4 — Ch is the effective local moment. The 
impurity-induced correction to the local susceptibility is thus 



Ax.=^[AXq«! 



q«q 



X 



^(A.! 



q«q 



Avq)], (10) 



and can be expressed as Axi ~ ^x/^ + Xi^ i-S-^ 
as a sum of the perturbed uniform susceptibility Ax = 

NACq=a^q=o/T = ^I'^/T, whereas 



X(cj) = -[^'^l+ojM{uj)-S\-h^, 



(5) 



where S = rjx ^- Again, the local fluctuation-dissipation re- 
lation 



duj cth—x-^{uj) = Cij, 



(6) 



is used to fix x- 

The next step in an inhomogeneous system is to diago- 
nalize the correlation matrix Cy = 12\'~^>^{'^x)*'^x where 
CvA — Cava. From the sum rule (|6jl it then follows that 
one can simultaneously diagonalize the susceptibility matrix 
X»i(w) = EaXa(^)K)*w^ with 



Xa(w) = --qx/ii^^ + ij\uj - Sx), 



q 



o«qri 



Cql^O 



(11) 



It is evident that Ax corresponds to a free spin with the mo- 
ment /i ^ 1/2 (at low doping) introduced by the impurity. 
Nontrivial is the spatial distribution of Xi- Since the main 
contribution in Eq. ( Ill> arises from q ^ Q and Cq ^ /i^, 
we obtain from Eq. ( Ill> Xi ~ ~Xov i-^-' ^^^ local response 
is just the intersite susceptibility of a homogeneous system, 
being the Fourier transform of Xq. 

In Fig. 3a we present results for 1/xi for the n.n. site. We 
use Eq. dlH with Xq determined via a self-consistent solution 
to Eqs. Q,(|2|i. To be consistent with the t-J model on a 2D 
square lattice we take here the form (7° = al{K^ + Cq) — b 
(7) with Cq = 2{cos{qx) + cos{qy) + 2) whereas a, b are chosen 



such that C°^g = and Cf^ == 1/2. As in the case of ho- 
mogeneous Xq in Fig. 1, we notice that the behavior is close 
to the Curie-Weiss form 1/xi oc T + 0i with a quahtative 
transition from 0i ~ for k < 0.7 to finite and large 8i 
for K > 1. Deviations from the Curie-Weiss dependence are 
understandable since Eq. dlU includes contributions from all 
q where 8q > 8q. Nevertheless the behavior at q ~ Q is 
dominant since considered k are quite large. It should be re- 
minded that the essence of the theory of anomalous scaling of 
spin response in underdoped cuprates 1 15] lies in the fact that 
anomalous uj,T dependence appears for T > Tk ^ even 
for substantial and T-independent k. 
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latter affect mainly the magnitude but not the T-dependence 

ofXi. 

In summary, we have presented a theory for the magnetic 
response of nonmagnetic impurities in doped AFM. It is based 
on the assumptions, and also evidence, that certain quantities 
are not substantially modified in the vicinity of an impurity. 
This is in particular the case of correlations C^ , but also of 
the collective mode amplitude rj and damping 7. Note that 
such assumptions would not be valid within a normal Fermi 
liquid but are rather the consequence of strong correlations, 
i.e., Cij exhibit a kind of charge-spin separation. 

Within the present theory the local spin response around the 
impurity clearly reflects the one in the homogeneous system 
which is also anomalous. Or alternatively, the measurements 
of the impurity-induced susceptibility allow for the recon- 
struction of the bulk Xq and corresponding correlation length 
^, as has been established in some other systems with the spin 
gap EH. 

We have also shown that xq in a doped AFM follows quite 
well the Curie- Weiss behavior Consequently also local Xi 
exhibit similar behavior, although with some deviations due 
to contributions from q 7^ Q. In any case, deviations from 
Xi ^ Ai/{T + <di) diminish as we consider further neigh- 
bors or when the perturbation by the impurity is less local (or 
weaker). From the theory it is clear that Ai is related to AFM 
correlations Cqj. Even more important, the transition from 
the regime with Kondo scale 8i ^ to a finite and fast in- 
creasing Qi reflects the crossover from a non-Fermi liquid to 
a more normal Fermi-liquid regime. 

Authors acknowledge the support of the Ministry of Edu- 
cation, Science and Sport of Slovenia under grant PI -0044. 



Figure 3: Local susceptibility 1/xi for the site n.n. to the impu- 
rity vs. T (both in units of i) for different k: a) evaluated via the 
perturbative calculation, b) with the diagonalization procedure. 

Alternatively, we can solve the impurity problem as defined 
by Eqs.(|5}-Q directly, without relying on the perturbation ex- 
pansion. We first find eigenvalues Cx and corresponding v\ 
numerically. For the homogeneous doped system we assume 
that C? correspond to 2D-lattice C" as discussed above. The 
final results for local Xi, as obtained from Eq. (|8} for the same 
parameters as in the perturbation calculation, are shown in 
Fig. 3b. Note that results presented in Figs. 3a,b are quite sim- 
ilar, even quantitatively. However, the full calculation yields 
results even closer to the Curie-Weiss behavior. 

Of interest is also the spatial distribution of the local sus- 
ceptibilities Xi around the impurity. It is evident from Eq. dlU 
that it reflects the bulk Xoi • Moreover, if we consider only the 
high-T local Curie constant Ai, it follows from Xq ^ C^/T 
that Ai ~ — C'oi- Thus Ai measures the equal-time AFM cor- 
relations, changing sign in accordance with the sign of Cg^. It 
is also plausible that the local response Xi is quite sensitive to 
any deviations of Cij — C° around the impurity. However, the 
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